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ABSTRACT 



(N 

Context. EK Eri is one of the most slowly rotating active giants known, and has been proposed to be the descendant of 
a strongly magnetic Ap star. 

Aims. We have performed a spectropolarimetric study of EK Eri over 4 photometric periods with the aim of inferring 
f») the topology of its magnetic field. 

Methods. We used the NARVAL spectropolarimeter at the Bernard Lyot telescope at the Pic du Midi Observatory, along 
with the least-squares deconvolution method, to extract high signal-to-noise ratio Stokes V profiles from a timeseries 
of 28 polarisation spectra. We have derived the surface-averaged longitudinal magnetic field Bg. We fit the Stokes V 
profiles with a model of the large-scale magnetic field and obtained Zeeman Doppler images of the surface magnetic 
, strength and geometry. 

Results. Be variations of up to about 80 G are observed without any reversal of its sign, and which are in phase 
with photometric ephemeris. The activity indicators are shown to vary smoothly on a timescale compatible with the 
rotational period inferred from photometry (308.8 d.), however large deviations can occur from one rotation to another. 
The surface magnetic field variations of EK Eri appear to be dominated by a strong magnetic spot (of negative polarity) 
which is phased with the dark (cool) photometric spot. Our modeling shows that the large-scale magnetic field of EK 
Eri is strongly poloidal. For a rotational axis inclination of i = 60°, we obtain a model that is almost purely dipolar. 
Conclusions. In the dipolar model, the strong magnetic/photometric spot corresponds to the negative pole of the dipole, 
which could be the remnant of that of an Ap star progenitor of EK Eri. Our observations and modeling conceptually 
support this hypothesis, suggesting an explanation of the outstanding magnetic properties of EK Eri as the result of 
f^*) ' interaction between deep convection and the remnant of an Ap star magnetic dipole. 
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1. Introduction the period may change from season to season, an average 

sjl ■ period of 308 ± 2.5 d provides an acceptable explanation 

<ti ' S„ S l^™ ting aCti 7 e ? Pf 1 * ? u ( , R of the data. DBSS10 report a photometric ephemeris based 

1362, HD 27536) is considered as a candidate for being the Qn thig periodj confirm that the main light var i ations are 

descendant of a strongly magnetic Ap star ( Stepien 1993, due t(J dark spotg and tentatively suggest that the star is a 
Strassmeier et aJ. 1999). The detection of a strong surface di p i e -dominated oblique rotator viewed close to equator- 
magnetic field, dominated by a poloidal component, is con- Qn Baged Qn tMg picture they propoge that thg rotation 
sistent with this hypothesis (Auriere et al 2008). EK En period ig twice the pnotometric per i d. 
has now been surveyed photometrically during more than 
30 years and its magnetic activity studied spectroscopically 

(Dall et al. 2010, hereafter DBSS10). These authors demon- Following the initial magnetic observations presented 
strate that, while the photometric variation and possibly b y Auriere et al. (2008), we have now monitored EK Eri 
throughout 4 seasons, corresponding to about 4 photomet- 
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* Based on data obtained using the Telescope Bernard Lyot we present our new data and the associated analysis, which 

at Observatoire du Pic du Midi, CNRS/INSU and Universite de enables us to establish the rotational period of EK Eri and 

Toulouse, France. to model its surface magnetic topology. 
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2. Observations 

2.1. Observations with NARVA L 

The 28 observations of EK Eri were obtained between 
20 September 2007 and 21 March 2011 at the 2m tele- 
scope Bernard Lyot (TBL) of Pic du Midi Observatory. We 
used the NARVAL spectropolarimeter, which is a copy of 
the ESPaDOnS instrument at the Canada-France-Hawaii 
Telescope (Donati et al. 2006a). NARVAL was used in po- 
larimetric mode with a spectral resolution of about 65000, 
in the same configuration and with the same procedure as 
described by Auriere et al. (2008). 

To complete the Zeeman analysis, least-squares decon- 
volution (LSD, Donati et al. 1997) was applied to all ob- 
servations. We used a line mask calculated for an effective 
temperature of 5250 K, logg = 3.5, and a microturbulence 
of 2.0 km s , consistent with physical parameters reported 
by DBSS10. In the present case this method enabled us to 
average about 11,000 spectral lines and to derive Stokes V 
profiles with a S/N that improved by a factor of about 40 
in comparison with that for single lines. We then computed 
the longitudinal magnetic field Be in Gauss, using the first- 
order moment method (Rees & Semel 1979, Donati et al. 
1997). 

Spectra obtained with NARVAL enable us to simultane- 
ously obtain magnetic data, measurements of activity indi- 
cators as well as radial velocity (RV) measurements, so that 
studying correlations between these quantities is straight- 
forward. The activity of the star was monitored by using 
line-activity indicators measured from the NARVAL spec- 
tra. We computed the S'-index (defined from the Mount 
Wilson survey, Duncan et al. 1991) for the chromospheric 
Ca II H & K line cores. Our procedure was first calibrated on 
the main sequence solar- type stars of Wright et al. (2004) , 
then we added 0.03 to the index to fit measurements of 5 
giant stars observed by Duncan et al. (1991) and Young 
et al. (1989). We also measured the relative intensity with 
respect to the continuum (R c ) of the Ca n IR triplet (854.2 
nm component) and Ha. 

We measured the RV of EK Eri from the LSD Stokes / 
profiles using a gaussian fit. 

Table 1 reports the date, HJD, photometric phase 
(based on the ephemeris of DBSS10) , the Ca II H & K 
S'-index, R c values for Ca II 854.2 nm and Ha, the mea- 
sured Be values and their uncertainties, and the inferred 
radial velocities. 



2.2. Results: variations of Be 

Table 1 and Fig. 1 show that Bg does not reverse its polarity 
during the 4 observed seasons. Figure 1 displays Bg mea- 
surements with respect to HJD. Vertical marks at the top of 
Fig. 1 indicate the times of photometric minima according 
to the ephemeris of DBSS10, i.e. when the cool spot is cen- 
tred on the visible hemisphere of the star (phase 0.0). We 
see that the unsigned Be is approximately at maximum at 
these dates. In the third season (2009-2010), which is well- 
centred on the spot passage, \Bg\ is confirmed to peak at 
phase 0. This trend is clearly visible in Fig. 2 (upper panel) 
which contains the same data as in Fig.l, but folded with 
respect to the photometric ephemeris (DBSS10, P — 308.8 
days, HJDo = 245 3372.679). This is consistent with a pic- 
ture in which the center of the dark spot of EK Eri coin- 
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Fig. 1. Variations of the longitudinal magnetic field Bg in G 
versus HJD (minus 2450000). Different symbols represent 
different observing seasons: stars: 2007-2008, crosses: 2008- 
2009, triangles: 2009-2010, diamonds: 2010-2011. Vertical 
marks at the top of the figure indicate the times of photo- 
metric minima according to the ephemeris of DBSS10. 



cides with the strongest \Bg\. DBSS10 suggested that the 
dark spot could correspond to the pole of the remnant of 
the dipole of an evolved Ap star. 

Since Bg does not reverse its polarity, in this scenario, 
only one pole of the putative dipole is observed, and the 
rotational period would then be equal to the photometric 
period, and not its double, contrary to the suggestion of 
DBSS10. 

Figure 2 also shows that the (year to year) scatter of 
the points is significant with respect to the error bars. It 
could be due to an incorrect period, or in changes in activity 
level, as observed in photometry. In the scenario where the 
rotational period is the photometric one, this would show 
that an oblique rigid rotator model (ORM, Stibbs 1950) 
alone could not explain the variations of Be- Some dynamo 
process would therefore need to be invoked. 

2.3. Results: variations of activity indicators and RV 

Figure 2 (lower panels) shows the variations of the activ- 
ity indicators versus photometric phase. Though the points 
are scattered, their lower envelope follows the variations of 
\Bg\ with respect to phase, with a maximum near phase 0, 
and a minimum near phase 0.5. This trend is more visible 
for Ca II (5-indcx and IR lines) than for Ha. Some points 
are significantly above the lower contours of the plots and 
they are discussed below. In particular, during the 2010- 
2011 season, the modulation of activity indicators appears 
to disappear (this season is highlighted in red in Fig. 2). 
Error bars are computed for each measurement of Ca II IR 
and Ha measurement. They are generally of order 0.001, 
and never reach 0.01. As to 5-index, an internal error of 
0.007 is estimated for the whole sample (see also Sect. 2.4) 
in agreement with the day-to-day scatter (Table 1, 12/13 
November 2007, 19/20 January 2008). These error bars for 
activity indicators are smaller than the symbols used in Fig. 
2 (lower panel) and are therefore not plotted there. 

Figure 3 shows the variations of RV versus phase. The 
long term stability of NARVAL is about 30 m s _1 (e.g. 
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Table 1. Log of observations of EK Eri (for details, see Sect. 2.1) 
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Date 




HJD 


Phot. 


S'-index 


Ca n 


Ha 


B ; 


a 


RV 








2450000+ 


Phase 




854.2 nm 




G 


G 


km s _1 


20 


Sep. 


07 


4364.69 


0.212 


0.501 


0.319 


0.226 


-98.6 


1.0 


7.118 


12 


Nov. 


07 


4417.44 


0.383 


0.457 


0.291 


0.233 


-21.1 


0.7 


7.047 


13 


Nov. 


07 


4418.45 


0.386 


0.456 


0.297 


0.233 


-21.4 


0.9 


7.053 


19 


Jan. 


08 


4485.31 


0.603 


0.470 


0.303 


0.230 


-28.9 


1.2 


7.096 


20 


Jan. 


08 


4486.38 


0.606 


0.463 


0.295 


0.227 


-31.6 


0.9 


7.110 


06 


Feb. 


08 


4503.35 


0.661 


0.469 


0.302 


0.234 


-45.6 


0.8 


7.111 


03 


Apr. 


08 


4560.31 


0.845 


0.513 


0.339 


0.250 


-62.7 


1.2 


7.145 


28 


Aug. 


08 


4707.68 


0.323 


0.502 


0.331 


0.251 


-12.8 


0.7 


7.052 


14 


Sep. 


08 


4724.60 


0.378 


0.511 


0.338 


0.262 


-27.9 


1.4 


7.043 


24 


Sep. 


08 


4734.57 


0.410 


0.510 


0.341 


0.270 


-18.8 


1.1 


7.000 


29 


Sep. 


08 


4739.71 


0.426 


0.506 


0.338 


0.264 


-13.7 


0.8 


6.996 


20 


Dec. 


08 


4821.37 


0.691 


0.543 


0.365 


0.266 


-36.8 


1.1 


7.046 


30 


Jan. 


09 


4862.35 


0.824 


0.588 


0.390 


0.291 


-46.4 


1.0 


7.207 


25 


Feb. 


09 


4888.31 


0.908 


0.573 


0.369 


0.270 


-67.8 


0.9 


7.176 


28 


Sep. 


09 


5103.62 


0.605 


0.602 


0.423 


0.286 


-48.0 


1.2 


7.080 


27 


Oct. 


09 


5132.57 


0.699 


0.576 


0.396 


0.280 


-56.9 


1.3 


7.101 


15 


Jan. 


10 


5212.44 


0.957 


0.547 


0.339 


0.254 


-88.6 


2.7 


7.190 


03 


Feb. 


10 


5231.32 


0.018 


0.566 


0.354 


0.259 


-90.7 


2.1 


7.200 


13 


Feb. 


10 


5241.29 


0.051 


0.575 


0.363 


0.263 


-90.9 


2.1 


7.214 


06 


Mar. 


10 


5262.30 


0.119 


0.548 


0.353 


0.258 


-77.0 


1.2 


7.218 


22 


Mar. 


10 


5278.32 


0.171 


0.533 


0.348 


0.260 


-58.9 


1.0 


7.195 


21 


Sep. 


10 


5461.62 


0.764 


0.598 


0.381 


0.278 


-90.9 


1.2 


7.062 


13 


Oct. 


10 


5483.64 


0.836 


0.565 


0.370 


0.262 


-94.6 


0.9 


7.107 


12 


Nov. 


10 


5513.55 


0.932 


0.571 


0.365 


0.264 


-93.1 


1.1 


7.183 


04 


Dec. 


10 


5535.36 


0.003 


0.576 


0.359 


0.270 


-87.2 


6.8 


7.150 


03 


Jan. 


11 


5565.28 


0.100 


0.585 


0.374 


0.268 


-77.7 


1.0 


7.150 


23 


Jan. 


11 


5585.44 


0.165 


0.589 


0.436 


0.272 


-57.6 


1.0 


7.166 


21 


Mar. 


11 


5642.31 


0.350 


0.574 


0.403 


0.277 


-39.5 


1.2 


7.130 



Note: Date of mid-observation, corresponding HJD, photometric phase, S'-index for Ca II H & K, R c values for Ca II 854.2 nm 
and Ha, Be values and their uncertainties, and radial velocity of the LSD Stokes / profile. 



Auriere et al. 2009). RV is found to vary with photometric 
phase, as already reported by Dall et al. (2005, DBSS10). 
There is some scatter near phase and a delay of the RV 
maximum by about 0.1 cycles with respect to Bg. This 
trend strongly suggests that the cool/magnetic spot is re- 
sponsible for the RV variations (see also Dall et al. 2005, 
DBSS10). Because of the small vsini of EK Eri, classi- 
cal Doppler imaging is not possible for this star. On the 
other hand, a small v sin i does not prevent mapping using 
Zeeman-Doppler imaging as described in Sect. 3. (e.g. Petit 
et al. 2008). 



2.4. Discrepant points in Bi and activity indicators 

Figure 2 shows that the seasonal variations of Bi and the 
activity indicators generally vary coherently with photo- 
metric phase, and that there is seasonal scatter greater than 
the measurement error bars. In addition there are some out- 
standing deviations from the mean variation. 

As far as the Bi measurements are concerned, the ob- 
servations of 20 September 2007 (the first one, which is our 
"discovery observation"), as well as those obtained on 21 
September and 13 October 2010 present outstanding devi- 
ations from the general phase variations. On 20 September 
2007, the S'-index and other activity indicators do not de- 
viate from the general phase variation. On the other hand, 
for the 2 first autumn 2010 observations, both Bi and activ- 
ity indicators are stronger than expected at their respective 
phases. 



Remarkably, in the 2010-2011 season, the modulation 
of activity indicators disappears: this season is highlighted 
in red in Fig. 2. While Bi values approach the typical level 
and modulation at the end of 2010, the intensity of the ac- 
tivity indicators remains high. Also, on some occasions, in 
particular 28 September 2009, the activity indicators are 
surprisingly high, while the Bi value is quite normal with 
respect to the phase. Looking closely at the data, we find 
that we may have a problem of normalization of the con- 
tinuum for this date, which can affect the measurements of 
activity indicators, but not Bg (as the Stokes V measure- 
ment is a differential measurement). To further investigate 
this potential background problem, we have computed the 
Ca II H emission index used by Morgenthaler et al. (2011, 
2012). This method makes use of a synthetic spectrum from 
the POLLUX database (Palacios et al. 2010) to normalise 
the continuum. It was found to be very effective in the case 
of solar- type dwarfs, enabling those authors to reach inter- 
nal errors of about 0.001 for chromospheric emission index 
measurements (Morgenthaler et al. 2011, 2012). However, 
in the case of EK Eri this method did not allow us to reduce 
the deviations, and we present our S'-index measurements 
in Table 1 and Fig. 2 as they are. 

Of peculiar interest are the isolated enhancements of 
Bg, activity indicators or both. They may be associated 
with flares as observed in active giants (e.g. Konstantinova- 
Antova et al. 2000, 2005). In this case it would be the first 
time that simultaneous observations of magnetic field and 
activity indicators have been performed during flares in the 
stellar context. In the solar simultaneous increase of 
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Fig. 2. Variations of the longitudinal magnetic field Bi (up- 
per panel) and of the activity indicators (lower panels) , with 
photometric phase. The lower panels, from top to bottom, 
show: S- index, R c for Ca 11 IR lines and H Q . Different sym- 
bols represent different observing seasons: stars: 2007-2008, 
crosses: 2008-2009, triangles: 2009-2010, diamonds: 2010- 
2011; the 2010-2011 season is shown in red. 
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Fig. 3. Variations of radial velocity with photometric 
phase. An error bar of 30 m is indicated. Different 
symbols represent different observing seasons: stars: 2007- 
2008, crosses: 2008-2009, triangles: 2009-2010, diamonds: 
2010-2011; the 2010-2011 season is shown in red. 



magnetic field and Ha has been observed (Lozitsky et al. 
2000). 

The outlying magnetic observations of September 2007 
and autumn 2010 were found to increase dramatically the 
reduced \ 2 01 our model fitting in Sect. 3. We will use only 
data between November 2007 and March 2010 (spanning 
about 3 rotations) for forthcoming magnetic analysis. 

3. The topology of the magnetic field at the 
surface of EK Eri 

3.1. Zeeman-Doppler imaging 

In order to model our series of Stokes V line profiles and to 
reconstruct the surface magnetic geometry of the star, we 
have used the Zeeman-Doppler imaging inversion method 
(ZDI; Donati and Brown 1997). The version of the code em- 
ployed is that described by Donati et al. (2006b): the sur- 
face magnetic field is projected onto a spherical harmonics 
frame and the magnetic field is resolved into poloidal and 
toroidal components. Because of the very slow rotation of 
EK Eri, its wsin? has been impossible to determine up to 
now due to the degeneracy between vsini and macroturbu- 
lence (DBSS10). These authors consider that 0.8 km s _1 is 
a safe upper limit for v sin i. In our models we used v sin i = 
0.5 km s _1 and found that varying this value did not signif- 
icantly change our results. We used a linear limb darkening 
coefficient equal to 0.75. We limited the spherical harmon- 
ics expansion to £ < 3 since increasing this threshold do 
not significantly change the results. 



3.2. The rotational period of EK Eri 

Knowing the rotational period of EK Eri is essential to be 
able to perform ZDI. Photometric variations of EK Eri have 
been now monitored during more than 30 years and a period 
of 308.8 ± 2.5 days has been adopted (DBSS10) as well as 
an ephemeris. This period may have changed slightly during 
the span of the photometric observations (Strassmeier et al. 
1999). DBSS10 argued that the rotational period could be 
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equal to twice the photometric period, i.e. 617.6 ± 5 days. 
To determine independently the rotational period of EK 
Eri, we have followed the approach of Petit et al. (2002), 
and calculated a set of magnetic maps, assuming for each 
map a different value for the rotational period. We impose 
a constant entropy for all the images and calculate as a 
goodness-of-fit parameter the reduced x 2 (Xr hereafter) by 
comparing the set of synthetic Stokes V profiles produced 
by ZDI to the observed time-series of Stokes V profiles. 
We then study the variations of x 2 over the range of rota- 
tion periods, to determine the period value producing the 
best magnetic model (identified by the lowest value of xV>- 
Here we scanned 400 values of the period between 200 and 
700 days, a range which encompasses the photometric and 
rotational periods proposed for EK Eri by DBSS10. Our pe- 
riodograms show a favoured period near 311 days. However 
the minimum x 2 value is high (about 25), notwithstanding 
that we limited our observational sample to the timespan 
November 2007 - March 2010 in order to avoid the effect of 
the strong deviations observed on 20 September 2007 and 
during the autumn of 2010. We therefore introduced a pos- 
sible differential rotation (Petit et al. 2002). Ultimately, we 
could not improve the results and deduced a negligible dif- 
ferential rotation for EK Eri. We also considered that our 
approximate line profile model might be the source of the 
high xl'- the current local line profile used in our model is 
a gaussian function. In an attempt to improve the fit of the 
line profile, we introduced a Lorentzian function and some 
asymmetry, but we still could not reduce significantly the 
X 2 - Ultimately, we consider that the relatively poor detailed 
fit to the Stokes V profiles resulting from our best-fit model 
is mainly due to the seasonal scatter of the magnetic field 
strength. 

Figure 4 shows the periodogram obtained for the whole 
data set between November 2007 and March 2010, assum- 
ing a rotational axis inclination of EK Eri of i = 60° and 
truncating the expansion at I = 3. For rotational periods of 
311.5 d, 308.8 d, and 625 d, we obtain respectively values for 
xl of 25 (our minimum), 25.3 and 38.2. The 600 day period 
suggested by DBSS10 is therefore eliminated both by our 
periodogram analysis and by the fact that their scenario to 
explain P rot = 2P p hot is not supported by the Be variations 
(Sect. 2.2). Since we found that phase 0.0 of the ephemeris 
of DBSS10 corresponds to approximately the center of the 
magnetic spot which dominates the magnetic map (see be- 
low), we decided to use the photometric period of 308.8 
days (and HJD = 245 3372.679) for our analysis of the 
magnetic topology. 

Figure 5 presents the fit of our best model to the ob- 
served Stokes V LSD profiles, for P Iot — 308.8 d, and 
i = 60°. 



3.3. Possible topologies of the surface magnetic field of EK 
Eri 

In the two previous subsections, we have presented the 
models and the parameters employed to infer the magnetic 
topology of EK Eri. Because vsmi cannot be determined 
precisely, the value of the inclination i is also only weakly 
constrained. This is our main remaining uncertain param- 
eter which we have attempted to constrain as part of the 
ZDI reconstruction. We have therefore fit our Stokes V data 
with our models and i varying from 85° to 40°. 




200 300 400 500 600 700 

Rot. Period (d) 



Fig. 4. Periodogram obtained using the entire data set be- 
tween November 2007 and March 2010. 



For all models we find that the poloidal component con- 
tains more than 98% of the reconstructed magnetic energy. 
When i decreases, we find the ZDI procedure to provide us 
with a better magnetic model, in the sense that the total in- 
formation content of the reconstructed magnetic geometry 
decreases (we follow here the maximum entropy criterion, 
Donati 2001, explicited in Morin et al. 2008). For high incli- 
nation, i > 80°, the topology is dominated by a quadrupo- 
lar configuration. In this case, the two rotation poles corre- 
spond to positive polarity while there is a nearly equatorial 
band, with a single enormous magnetic spot, both of nega- 
tive polarity. The magnetic spot corresponds to the phase 
of the dark (and therefore presumably cool) photometric 
spot to better than 0.1 cycles. 

For i < 80°, the dipole component dominates and cor- 
responds to about 90% of the magnetic energy at about i 
= 60°. Figure 6 shows the corresponding ZDI map. This 
map also features a strong magnetic spot of negative po- 
larity, corresponding to the negative pole of the magnetic 
dipole. Again the magnetic maximum corresponds in phase 
to the photometric spot. The characteristics of this model 
are compatible with the scenario of Auriere et al. (2008) and 
DBSS10, in which the magnetic/photometric spot would 
correspond to the remnant of a magnetic pole of the Ap 
star progenitor of EK Eri. 

Given the lower information content (i.e. higher en- 
tropy) associated with the predominantly dipolar config- 
uration, we consider this model to be better, or at least 
more likely, than the quadrupolar model. In addition, this 
model suggests an obvious physical interpretation: that the 
corresponding magnetic pole is at the position of the pho- 
tometric spot of EK Eri, and that the dipole could be the 
remnant of the magnetic field of an Ap star progenitor. 
Both Fig. 5 - presenting the variations of the Stokes V pro- 
files - and Fig. 2 (upper panel) - presenting the variations 
of Be - show that the magnetic field varies on a large scale, 
and does not experience a sudden drop when the "mag- 
netic spot" is not visible, as it would do in the case of 
a very localized strong magnetic spot: a dipolar configu- 
ration naturally explains the observed behavior. However, 
even if the quadrupolar configuration is less favoured by 
our modeling procedure and if this topology does not ap- 
peal for an obvious formation scenario, we cannot rule it 
out as a possible solution. As explained in Sect. 2.4, the 
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Fig. 5. Fit of our best-fit model (for P rot = 308.8 d, i = 
60°) (thin red curves) to the Stokes V LSD profiles (thick 
black curves), after correction of the mean radial velocity 
of the star. Rotational phases (according to the ephemeris 
of DBSS10) are indicated. 



present ZDI study was made using only the 20 spectra cor- 
responding to about 3 rotations, between November 2007 
and March 2010, to avoid epochs when activity appeared 
to deviate significantly form the mean variation. However, 
the main results as described above would be the same if 
we employed the entire data set, except that the Xr would 
have been significantly worse. 

From the dipolar model (corresponding to i = 60° 
and Fig. 6) we can measure /3 (the angle between rota- 
tional and magnetic axes of the dipole component) and the 
dipole strength (since the dipolar component dominates, 
the strength of the dipole corresponds to about the max- 
imum magnetic field measured on our map): we find j3 — 
45°, B d ~ 207 G. 



We have also fit our data with a model with a pure 
dipole, corresponding to spherical harmonic I = 1 (and 
assuming i = 60°). The goodness of fit is (naturally) worse 
than for the models described above, and the inferred dipole 
parameters are /3 = 22°, B d ~ 198 G. 

The magnetic maximum is measured to occur at phases 
0.04 and 0.96, respectively, for the two ZDI models de- 
scribed above. 

It could also be of interest to use the measurements of 
Bi from Sect. 2.2 (Table 1) to fit a dipolar oblique rotator 
model using the relations of Preston (1967) and reviewed by 
Auriere et al. (2007). In this case we find (3 = 23° and dipole 
strength of 370 G (taking i = 60° as suggested above). 

The smaller value obtained for B d from our ZDI models 
with respect to B d inferred from the Bi variations is due in 
part to the fact that our ZDI models correspond to Stokes 
V signatures at phases near magnetic maximum that are 
weaker than those observed (Figure 5). 

Taking i — 60°, in the case of rigid rotation, using R 
= 4.68 R© from the evolutionary model of Charbonnel & 
Lagarde (2010) used by Auriere et al. (2008), we get wsini 
= 0.69 km s -1 . 



4. Summary and discussion 

4.1. Summary 

Spectropolarimetric observations of EK Eri have been ob- 
tained during 4 seasons, corresponding to about 4 photo- 
metric periods of 308.8 d. duration. 

The longitudinal magnetic field Be is detected consis- 
tently, with Be variations of up to about 80 G without any 
reversal of its sign. The maximum of \Bg\ is observed to 
coincide with the phase of photometric light minimum (i.e. 
phase 0.0 of the photometric ephemeris of DBSS10). 

The activity indicators measured from Ca n H & K, the 
Ca IR triplet and Ha are shown to vary with small ampli- 
tude and in phase with \Bg\. However some strong devia- 
tions are observed, which may suggest episodes of flares. 

The radial velocity also varies in a reasonably coherent 
way according to the photometric ephemeris, but with a 
phase delay of about 0.1 cycles. 

For our Zeeman Doppler imaging (ZDI) investigation 
we used the 308.8 d photometric period (DBSS10) which is 
similar to our best rotational period, based on quality-of-fit 
of the ZDI maps, of about 311 d. 

For all the described models we find that the large scale 
magnetic field of EK Eri is poloidal (more than 98 % of 
the magnetic energy is contained in this component). In 
addition, the surface magnetic field variations of EK Eri 
appear to be dominated by a strong magnetic spot (of neg- 
ative polarity) which is phased with the darker, redder (and 
therefore presumably cool) photometric spot. 

The ratio between dipolar and quadrupolar components 
of the magnetic field increases with decreasing value of the 
rotational axis inclination i. Whereas the quadrupolar com- 
ponent dominates for i greater than 80°, for i = 60°, we 
obtain a model that is almost purely dipolar. In the dipo- 
lar model, the photometric spot at phase corresponds to 
the pole of negative polarity of the dipole, which could be 
the remnant of that of the Ap star progenitor of EK Eri. 
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Fig. 6. Magnetic map of EK Eri for i — 60°. Each chart 
illustrates the field projection onto one axis of the spherical 
coordinate frame. The magnetic field strength is expressed 
in G. The star is shown in flattened polar projection down 
to latitudes of —30° , with the equator depicted as a bold 
circle and parallels as dashed circles. Radial ticks around 
each chart indicate the observed rotational phases. 



4.2. Discussion 

The evolution of a magnetic (2 M Q ) Ap star from the main 
sequence to the red giant branch is characterized by a sur- 
face magnetic field which is expected to weaken (as 1/R 2 if 
conservation of magnetic flux is assumed) while the convec- 
tive envelope deepens. The rotational period of the surface 
also lengthens as the star expands. For an Ap star progeni- 
tor of EK Eri, one can infer a magnetic field of several thou- 
sand gauss and a very thin convective envelope on the main 
sequence (Stepien 1993, Auriere et al. 2008). At the main 
sequence evolutionary stage, the magnetic field would be 
sufficient to suppress convection (e.g. discussion by Theado 
et al. 2005). At the present evolutionary stage of EK Eri 
we measure a large scale surface magnetic field dominated 
(for i = 60°) by a dipole of about 200 G strength, when 
the convective envelope of our 2 Mq model contains ~ 0.37 
Mq of the star's mass (Auriere et al. 2008). The outstand- 
ing magnetic activity properties of EK Eri are therefore the 
result of the interplay of the remaining magnetic field from 
one Ap star and deep convection. Such an interaction be- 
tween a pre-existing magnetic field and thermal convection 
has been studied through numerical simulations mostly in 
the solar context (see e.g. Hurlburt et al. 1996, Cattaneo 
et al. 2003, Proctor 2004, Thomson 2005, Strugarek et al. 
2011) and in the cases of core convection (Featherstone et 
al. 2009). For example, the transition between regimes when 
the imposed magnetic field is strong (magnetoconvection) 
and when the magnetic field is weak (dynamo) is investi- 
gated (e.g. Cattaneo et al. 2003). 

No differential rotation is inferred from our modeling 
(Sect. 3.2): this can be expected both from a strong (fossil) 
magnetic field and from a slow rotator. Figure 5 shows that 
an average magnetic field can be defined during the 3 ob- 
served rotational periods used for modeling. Figure 2 shows 
that significant deviations occur during that period of time, 
and some are even stronger before (20 September 2007) 
and afterward (21 September and 13 October 2010; see also 
Sect. 2.4). This shows that even if a magnetic dipole domi- 
nates the magnetic topology, we are not in the presence of 
a stable, large scale magnetic field as in an Ap star. Figure 
5 also shows that it was not possible to perfectly fit our 
observations at phases near maximum magnetic strength. 
Therefore, the interplay of the remaining field and convec- 
tion appears to induce some variable magnetic component. 

Another striking difference between the magnetic spot 
of EK Eri and the magnetic pole of an Ap star, is that for 
an Ap star surface features ("spots") correspond generally 
to chemical over-abundances of some elements, and are ex- 
pected to be bright (Krticka et al. 2007, Liiftinger et al. 
2010). In the case of EK Eri, the spot is dark and redder 
(and presumably cool), while the abundances are found to 
be similar to those of the sun at all phases (DBSS10; Van 
Eck et al. 2011 in preparation); the spot of EK Eri there- 
fore appears to be a temperature spot, as observed in the 
late-type stars, where dynamo operates. As discussed by 
Auriere et al. (2008) and Van Eck et al. (2011, in prepa- 
ration), the overabundances observed at the surface of an 
Ap star would have been erased by the deep mixing occur- 
ing during the first dredge-up in which EK Eri is currently 
engaged. The differences between models for the magnetic 
poles of Ap stars and sunspots have been investigated (e.g. 
Shibahashi 2004). 
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EK Eri therefore appears to be an excellent laboratory 
to study the interplay between magnetic field and convec- 
tion: the magnetic field at the surface of EK Eri may be 
considered as a new type of magnetic field, neither a rigid 
large scale magnetic field (e.g. dipole and oblique rotator 
model) nor a solar-type dynamo magnetic field. This mag- 
netic field has the potential to strongly interact with the 
oscillations detected in EK Eri (Dall et al. DBSS10, 2011). 
However, the dipole field measured with our modeling from 
2007-2010 data is rather weaker than the 300-500 G ex- 
pected from preliminary computations (Dall et al. 2011) to 
be able to significantly alter the pulsations. 

This work shows that the magnetic variations of EK Eri 
arc dominated by stellar rotation (i.e. rotational modula- 
tion of a dipole in our prefered scenario). The survival of 
a fossil dipole can explain the long term stability of the 
location of the dark spot, and magnetic cycles of varying 
photometric amplitude may be due to dynamo-like activ- 
ity. Because of the very small v sin i of EK Eri, only large 
scale features of the surface magnetic field can be mapped 
with ZDI. However, the activity indicators are sensitive 
to all scales of magnetic regions. The observations of the 
2010-2011 season show (Fig. 2) that while in September- 
October 2010 both the magnetic field and activity indica- 
tors were stronger than expected for that rotational phase, 
the magnetic field became "normal" at the beginning of 
2011 while the values of the activity indicators remained 
high. Therefore the large scale magnetic component could 
be more stable than the small-scale field . In this context, it 
would be worthwhile to investigate if the possible magnetic 
cycles suggested by photometry (DBSS10) exist as well in 
the magnetic data. Studying the correlation between photo- 
metric amplitude and surface magnetic strength/topology 
of EK Eri would certainly help to understand the interplay 
of the dipole remnant and convection. 
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